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We present the first measurement of dijet angular distributions in pp collisions at ^/s = 1.96 TeV 
at the Fermilab Tevatron Collider. The measurement is based on a dataset corresponding to an 

integrated luminosity of 0.7 fb~^ collected with the DO detector. Dijet angular distributions have 
been measured over a range of dijet masses, from 0.25 TeV to above 1.1 TeV. The data arc in good 
agreement with the predictions of perturbative QCD and are used to constrain new physics models 
including quark compositeness, large extra dimensions, and TeV~^ scale extra dimensions. For all 
models considered, we set the most stringent direct limits to date. 

PACS numbers: 12.60.Rc, 11.25. Wx, 12.38.Qk, 13.87.Ce 



At large momentum transfers, dijet production has the It can be used to test the standard model (SM) at previ- 
largest cross section of all processes at a hadron collider ously unexplored small distance scales and to search for 
and therefore directly probes the highest energy regime. signals predicted by new physics models. The angular 
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distribution of dijets with respect to the hadron beam 
direction is directly sensitive to the dynamics of the un- 
derlying reaction. While in quantum chromodynamics 
(QCD) this distribution shows small but noticeable de- 
viations from Rutherford scattering, an excess at large 
angles from the beam axis would be a sign of new physics 
processes not included in the SMjSuch as substructure of 
quarks ( "quark compositeness" ) [H, 0, 0] , or the existence 
of additional compactified spatial dimensions ( "extra di- 
mensions") [1, H, S 0, Si- Earher measurements of dijet 
angular distributions and related observables in pp colli- 
sions at -y/fi = 1.8 TeV were used to set limits on quark 
compositeness (il. [loj. 

In this Letter we present the first measurement of 
dijet angular distributions in pp collisions at a center- 
of-mass energy of i/i = 1.96 TeV. The data sample, 
collected with the DO detector during 2004-2005 in 
Run II of the Fermilab Tevatron Collider, corresponds 
to an integrated luminosity of 0.7 fb~^. In the experi- 
ment and in theory calculations, jets are defined by the 
Run II midpoint cone jet algorithm ll| with a cone 
radius oi TZ = (A?/)^ -I- (Ai/))^ = 0.7 in rapidity y 
and azimuthal angle (j). Rapidity is related to the po- 
lar scattering angle 6 with respect to the beam axis by 
y = O.51n[(l-h/3cos6')/(l-/3cos6i)] with /3 = \p\/E. 
We measure distributions in the dijet variable Xdijct = 
exp(|yi — J/2I) in ten regions of dijet invariant mass Mjj, 
where yi and 2/2 are the rapidities of the two jets with 
highest transverse momentum px with respect to the 
beam axis in an event. For massless 2 — > 2 scatter- 
ing, the variable Xdijct is related to the polar scatter- 
ing angle 6* in the partonic center-of-mass frame by 
Xdijet = (1 + cos6'*)/(l - cos^*). The choice of this 
variable is motivated by the fact that Rutherford scat- 
tering is independent of Xdijct- The phase space of this 
analysis is defined by Mjj > 0.25 TeV, Xdijet < 16, and 
2/boost = 0.5 \yi + 2/2I < 1- Together, the Xdijet and yboost 
requirements restrict the jet phase space to \yjot\ < 2.4 
where jets are well-reconstructed in the DO detector and 
the energy calibration is known to high precision. To 
minimize sensitivity to correlated experimental and the- 
oretical uncertainties, the Xdijct distributions in the dif- 
ferent Mjj ranges are normalized by their respective inte- 
grals. Based on the measurement, we set limits on quark 
compositeness [l|, 0, 01, large spatial extra dimensions 
according to the model proposed by Arkani-Hamed, Di- 
mopoulos and Dvali (ADD LED) Mlsl, and TeV~^ scale 
extra dimensions (TeV~^ ED) [y, lZ|, l8(. 

A detailed description of the DO detector can be found 
in Ref. [13]. The event selection, jet reconstruction, jet 
energy and momentum correction in this measurement 
follow closely those used in our recent measurement of 
the inclusive jet cross section [l3| . The primary tool for 
jet detection is the finely segmented uranium-liquid ar- 
gon calorimeter that has almost complete solid angular 
coverage 1.7° <0< 178.3° [l3|- Events are triggered by 



the jet with highest pT, referred to as p^^^. In each Mjj 
region, events are taken from a single trigger which is cho- 
sen such that the smallest p^'^^ in the Mjj region is above 
the threshold that ensures 100% efficiency. The My re- 
gions utilize triggers with different prescales, resulting 
in integrated luminosities of O.lOpb"^ {Mjj < 0.4 TeV), 
1.54 pb-i (0.4 < Mjj < 0.5 TeV), 17pb-i (0.5 < Mjj < 
0.6TeV), 73pb-i (0.6 < Mjj < 0.8 TeV), 0.5 fb-^ (0.8 < 
Mjj < 1.0 TeV), and 0.7fb-i {Mjj > 1.0 TeV). 

The position of the pp interaction is reconstructed 
using a tracking system consisting of silicon microstrip 
detectors and scintillating fibers, located inside a 2 T 
solenoidal magnet 12] , and is required to be within 50 cm 
of the detector center along the beam direction. The 
jet four-momenta are corrected for the response of the 
calorimeter, the net energy flow through the jet cone, 
energy from event pile- up and multiple pp interactions, 
and for systematic shifts in ]?/] due to detector effects [13]. 
Cosmic ray backgrounds are suppressed by requirements 
on the missing transverse momentum in an event [l3| . 
Requirements on characteristics of shower shape are used 
to suppress the remaining background due to electrons, 
photons, and detector noise that mimic jets. The effi- 
ciency for these requirements is above 97.5%, and the 
fraction of background events is below 0.1% in all Mjj 
regions. 

The Xdijct distributions are corrected for instrumental 
effects using events generated with pythia v6.419 (i3] 
using tune QW jH and MSTW2008LO parton distribu- 
tion functions (PDFs) [IB]- The generated particle-level 
events are subjected to a fast simulation of the DO de- 
tector response, based on parametrizations of resolution 
effects in px, the polar and azimuthal angles of jets, jet 
reconstruction efficiencies, and misidentification of the 
event vertex. These parametrizations have been deter- 
mined either from data or from a detailed simulation of 
the DO detector using GEANT [i3] . The generated events 
are reweighted according to the Mjj distribution in data. 
To minimize migrations between Mjj regions due to reso- 
lution effects, we use the simulation to obtain a rescaling 
function in Mjj that optimizes the correlation between 
the reconstructed and true values. The bin sizes in the 
Xdijct distributions are chosen to be much larger than the 
Xdijct resolution. The bin purity after Mjj rescaling, de- 
fined as the fraction of all reconstructed events that were 
generated in the same bin, is between 42% and 68%. We 
then use the simulation to determine Xdijet bin correction 
factors for the differential cross sections in the different 
Mjj regions. These also include corrections for the ener- 
gies of unreconstructed muons and neutrinos inside the 
jets. The total correction factors for the differential cross 
sections are typically between 0.9 and 1.0, and always in 
the range 0.7 to 1.1. The corrected differential cross sec- 
tions within each Mjj range are subsequently normalized 
to their integrals, providing the corrected, final results 
for 1 /(Xdijct • da /dxdijct at the "particle level" as defined 
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FIG. 1: Normalized differential cross sections in Xdijct com- 
pared to standard model predictions and to the predictions 
of various new physics models. The error bars display the 
quadratic sum of statistical and systematic uncertainties. The 
standard model theory band includes uncertainties from scale 
variations and PDF uncertainties (see text for details). 



in Ref. {1$. 

In order to take into account correlations between sys- 
tematic uncertainties, the experimental systematic un- 
certainties are separated into independent sources, for 
each of which the effects are fully correlated between all 
data points. In total we have identified 76 independent 
sources, of which 48 are related to the jet energy calibra- 
tion and 15 to the jet px resolution uncertainty. These 
are the dominant sources of uncertainty. Smaller con- 
tributions are from the jet 9 resolution and from the 
systematic shifts in y. All other sources are negligible. 
All sources and their effects are documented in Ref. [l^ . 
For Mjj < 1 TeV (Mjj > 1 TeV) systematic uncertainties 
are l%-5% (3%-ll%); they are in all cases less than the 
statistical uncertainties. 

The results are available in Ref. [l9| and displayed 
in Fig. [T] The normalized Xdijet distributions are 
presented in ten Mjj regions, starting from Mjj > 
0.25 TeV, and including one region for Mjj > 1.1 TeV. 
The data are compared to predictions from a pertur- 
bative QCD calculation at next-to-leading order (NLO) 
with non-perturbative corrections applied. The non- 
perturbative corrections are determined using PYTHIA. 



They are defined as the product of the corrections due 
to hadronization and to the underlying event. The 
NLO results are cornputed using FASTNLO [2^ based 
on NLO JET -I- -I- 21I, [221. All theory calculations use 
MSTW2008NLO PDFs [11] and the corresponding value 
of as{Mz) — 0.120. The PDF uncertainties are provided 
by the twenty MSTW2008NLO 90% C.L. eigenvectors. 
Renormalization and factorization scales fi are varied si- 
multaneously around the central value of /iq = (pt) in 
the range 0.5 /io < fJ- < 2 where (pt) is the average 
dijet Pt- The quadratic sum of scale and PDF uncer- 
tainties is displayed as a band around the central SM 
value in Fig. [TJ The scale (PDF) uncertainties are al- 
ways below 5% (2%) so the band is nearly a line. The 
theory, including the perturbative results and the non- 
perturbative corrections, is in good agreement with the 
data over the whole AIjj range with a (defined later) of 
127.2 for 120 data points in ten normalized distributions. 
Based on the agreement of the Xdijot measurement with 
the SM, we proceed to set limits on quark compositeness, 
ADD LED, and TeV'^ ED models. 

Hypothetically, quarks could be made of other parti- 
cles, as assumed in the quark compositeness model in 
Ref. [H, [3, We investigate the model in which all 
quarks are considered to be composite. The parameters 
in this model are the energy scale A and the sign of the 
interference term rj between the standard model and the 
new physics terms. The ADD LED model assumes 
that extra spatial dimensions exist in which gravity is al- 
lowed to propagate. Jet cross sections receive additional 
contributions from virtual exchange of Kaluza-Klein ex- 
citations of the graviton. There are two different for- 
mahsms (GRW and HLZ Q). The model parameter 
is the effective Planck scale, Ms, and the HLZ formalism 
also includes the subleading dependence on the number 
n of extra dimensions. The TeV"^ ED model 0, 0, S] 
assumes that extra dimensions exist at the TeV~^ scale. 
SM production cross sections are modified due to virtual 
Kaluza-Klein excitations of the SM gauge bosons. In this 
model, gluons can travel through the extra dimensions, 
which changes the dijet cross section. The parameter in 
this model is the compactification scale, Mc- 

The new physics contributions have only been calcu- 
lated to leading order (LO), while the QCD predictions 
are known to NLO. In this analysis, to obtain the best 
estimate for new physics processes, we multiply the new 
physics LO calculations bin-by-bin by the SM A:-factors 
{k ~ cTNLo/fLo)- The A:-factors are in the range 1.25- 
1.5, increasing with Aljj and decreasing with Xdijct- Their 
effects on single bins of the normalized Xdijet distribu- 
tions within the different Mjj regions is below 12%. The 
new physics cross sections are computed using the ma- 
trix elements from Refs. 0, 0, [H, l8!|. The theoretical 
variations (scale variations and PDF uncertainties) are 
consistently propagated into both the SM and the new 
physics contributions. Predictions for the different mod- 
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TABLE I: Expected and observed 95% C.L. limits in units of TcV on various new physics (NP) models for different Bayesian 
priors, and for the Ax^ criterion. 





Prior flat in 


NP Lagrangian 


Prior flat in 


NP x-section 




3.84 criterion 


Model (parameter) 


Expected 


Observed 


Expexted 


Observed 


Expected 


Observed 


Quark compositeness (A) 














?? = +1 


2.91 


3.06 


2.76 


2.84 


2.80 


2.92 




2.97 


3.06 


2.75 


2.82 


2.82 


2.96 


TeV-i ED (Mc) 


1.73 


1.67 


1.60 


1.55 


1.66 


1.59 


ADD LED (Ms) 














GRW 


1.53 


1.67 


1.47 


1.59 


1.49 


1.66 


HLZ n = 3 


1.81 


1.98 


1.75 


1.89 


1.77 


1.97 


HLZ n = 4 


1.53 


1.67 


1.47 


1.59 


1.49 


1.66 


HLZ n = 5 


1.38 


1.51 


1.33 


1.43 


1.35 


1.50 


HLZ n = 6 


1.28 


1.40 


1.24 


1.34 


1.25 


1.39 


HLZ n = 7 


1.21 


1.33 


1.17 


1.26 


1.19 


1.32 



els are compared to the Xdijet data and to the SM results 
in Fig. [1] It is observed that all models predict increased 
contributions as Xdijct 1 towards large Mjj. The Mjj 
evolution of the excess tovifards small Xdijct is observed 
to be different for different models. 

We define the b etween data and theory using the 
Hessian approach [25| which introduces nuisance param- 
eters for all correlated sources of experimental and the- 
oretical uncertainty. The x^ is then minimized with re- 
spect to all nuisance parameters, and is therefore only 
a function of the new physics model parameter (s). In 
most cases x^ ha-s the minimum for a new physics mass 
scale of infinity, corresponding to the SM value. Only for 
the quark compositeness model with positive interference 
and for the TeV^^ ED model x^ has small minima at 
A = 9.88 TeV with Ax^ = 0.01 and Mc = 2.96 TeV with 
Ax^ = 0.28 below the SM value, respectively. 

The x^ is then transformed into a likelihood which is 
used in a Bayesian procedure 10] to obtain 95% C.L. 
limits on the new physics mass scales A, Mc, and Ms in 
the different models. The prior is chosen to be flat in the 
new physics mass scale raised to the power in which it 
appears in the Lagrangian or, alternatively, raised to the 
power in which it enters the model cross section. While 
the former has been used in many previous analyses, the 
latter is statistically preferred for being unbiased in the 
cross section. Alternatively, we have applied a procedure 
which defines the 95% C.L. limit as the mass scale at 
which x^ — Xmin = 3-84 [2^. This procedure has the ad- 
vantage of being independent of an assumed prior. The 
observed limits and the expectation values are listed in 
Table m All observed limits are within one standard de- 
viation of the expected limits. 

The limit on Mc obtained in this analysis, while in- 
ferior to indirect limits from electroweak precision mea- 
surements (Ref. Q and references therein), is comple- 
mentary and is the result of the first direct search for 
TeV~^ extra dimensions at a particle collider. The lim- 
its on Ms in the different formalisms of the ADD LED 



model are on average slightly higher as compared to re- 
cent DO results from the combination of 1 fb~^ of dielec- 
tron and diphoton data in Ref. , which were so far the 
most restrictive limits on ADD LED. Our limits on quark 
compositeness improve previous results from related di- 
jet observables and are the most stringent limits 

to date. 

In summary, we have presented the first measurement 
of dijet angular distributions in Run II of the Fermilab 
Tevatron Collider. This is the first measurement of an- 
gular distributions of a hard partonic scattering process 
at energies above 1 TeV in collider-based high energy 
physics. The normalized Xdijct distributions are well- 
described by theory calculations in next-to-leading order 
in the strong coupling constant and are used to set lim- 
its on quark compositeness, ADD large extra dimensions, 
and TeV~^ extra dimensions models. For the TeV~^ ex- 
tra dimensions model this is the first direct search at a 
collider. For all models considered, this analysis sets the 
most stringent direct limits to date. 
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